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Abstract: In this article we discuss a set of geometric ideas which shed some 
light on the question of directed polymer pinning in the presence of bulk disor- 
der. Differing from standard methods and techniques, we transform the prob- 
lem to a particular dependent pcrcolative system and relate the pinning tran- 
sition to a percolation transition. 



1. Introduction 

Motivating example. The totally asymmetric exclusion process (TASEP for 
short) is defined as follows on the set Z: At time 0, a (possibly random) configuration 
of particles is given, in such a way that each site contains at most one particle. To 
each edge of the lattice is associated a Poisson clock of intensity 1. Whenever this 
clock rings, and there is a particle at the left-end vertex of this edge and no particle 
at the right-end vertex, the particle moves to the right; otherwise the ring of the 
clock is ignored. The product over Z of Bernoulli measures of density p G (0, 1) is 
invariant by this dynamics; in that case, the average number of particles passing 
through the origin up to time t is equal to p(l — p)t, i.e. the flux through a given 
bond is exactly p(\ — p). 

The process is modified at the origin, by imposing that the Poisson clock asso- 
ciated with the bond eo = (0, 1) is A > 0. When A > 1, one can still prove that 
the above expression for the flux holds asymptotically, although the Bernoulli mea- 
sure is not an equilibrium measure anymore. From now on, we shall assume that 
< A < 1. 

One of the fundamental questions in driven flow is to understand under which 
conditions such a static obstruction results in the formation of a "platoon" start- 
ing at the origin and propagating leftward. A convenient quantitative criterion for 
platoon formation is to start the TASEP with step initial conditions, i.e. all sites 
a; < filled and all sites x > 1 empty, and to consider the average current, j(A), 
in the long time limit t — > oo. j(0) = 0, j is non-decreasing, and j(A) = 1/4 for 
A = 1. Thus the issue is to determine the critical intensity Ac, which is defined as 
the supremum of all the A for which j(A) < j. Estimates for the value of Ac were 
given in |lOj and , and the full hydrodynamical picture is proved in [isj . 
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The blockage problem for the TASEP has been studied numerically and by exact 
enumerations. On the basis of these data, in the value Ac = 1 is conjectured. 
Recently this result has been challenged and Ac = 0.8 is claimed. 

One-dimensional driven lattice gases belong to the universality class of Kardar- 
Parisi-Zhang (KPZ) type growth models. In particular the asymmetric simple exclu- 
sion process can be represented as the so-called body-centered solid-on-solid version 
of (l-l-l)-dimensional polynuclear growth model, or as directed polymer subject to 
a random potential. The TASEP also has a well known representation in terms 
of last-passage percolation (or maximal increasing subsequences, known as Ulam's 
problem) . In this article wc choose to work in the setup of last passage percolation. 
The slow bond induces an extra line of defects relative to the disordered bulk. If 
A < Ac, the optimal path, i.e. the geodesic^ is pinned to the line of defects. As 
A — > Ac, the geodesic wanders further and further away from the line of defects and 
the density of intersections with the line of defects tends to zero. For A > Ac, the 
fluctuations of the geodesic are determined by the bulk, and the line of defects is 
irrelevant. 

In the present work we will not establish the actual value of Ac (since we cannot) 
or settle the question as to whether it is equal to 1; the main goal of this paper is to 
describe a new way of looking at the problem which gives some insight about the 
precise behavior of the system. More precisely, we show how the problem can be 
studied through particular dependent percolative systems constructed in such a way 
that the pinning transition can be understood in terms of a percolation transition. 



2. Interpretation as pinning in Ulam's problem 

In this section we describe the representation of our initial problem in terms of 




2.1. The Model 



Let p(^' be the distribution of a Poisson point process of intensity A*^^) = 1 in the 
plane R^, and let fi^^^ the set of all its possible configurations; for all n > 0, let 

('21 f2) 

Pn be the law of its restriction to the square Q„ := [0, n] x [0, n] and ri„ ' be the 
configuration space of the restricted process. 

The maximal increasing subsequence problem, or Ulam's problem, can be for- 
mulated in the following geometric way: Given a configuration cj^^' G fj'^^ and 
its restriction ui'if' to Q„, look for an oriented path tt (moving only upward and 
rightward) from (0,0) to (n,n) collecting as many points from ojn as possible; as 
in the case of last-passage percolation described in the introduction, we shall call 
such a path a geodesic. Let A/'„ denote the number of collected points along such 
an optimal path (which need not be unique). It is a well known fact (see 0) that 

lim t^J^ = 2. (2.1) 

n— »-|-oo n 

In a remarkable paper j^, K. Johansson showed that for this model the transver- 
sal fluctuations of the geodesies are of order n^^^. A closely related problem is 
considered in 0. 

We now modify the original model in the following way: Assume in addition, 
that on the main diagonal y ~ x there is an independent one-dimensional Poisson 
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point process of intensity A'-^-' = A, and let P^^-* be its distribution. We will denote 
by fl^^'' (resp. fil^'') the configuration space of this process (resp. of its restriction to 
Q„). Finally if w^^) is a realization of F^^^ and w^^) one of P^'^\ let lo = tjf^) y ^^(i) 
be their union, and LOn the restriction of lo to Q„. 

What can now be said about an optimal directed path starting at (0, 0) and 
ending at If A » 1, clearly the geodesic will stay close to the diagonal and 

lim -^^'^ =: e(A) > 2. (2.2) 

Instead of H2.2|l in our context it will be more convenient to use a more geometric 
notion. We say that the directed polymer is pinned (with respect to the diagonal) 
if p(^^ X P'^^^-almost surely, the number of visits made by geodesies to {(m, u) : < 
u < t} is of order jt, for some 7 > 0, for all t large enough. One expects the 
"energy" notion H2.2|l and geometric notion of pinning to be identical, but this is 
yet another point which remains to be proved. 

e(A) is non-decreasing and there is a critical value Ac, where e hits the value 2. 
The same arguments which predict Ac = 1 in the case of the TASEP yield Ac = in 
the case of our model. Thus any extra Poisson points along the diagonal are expected 
to pin the directed polymer. Such a behavior is extremely delicate, and the answer 
depends on the nature and behavior of the geodesies in the initial, unperturbed 
system. Very little is known, even on a heuristic level, when the underlying measure 
governing the behavior of the polymer is not "nice" (with a kind of Markov property, 
for example a simple symmetric random walk). Our criteria presented below give 
partial but rigorous answers as to whether Ac is strictly positive or not. 

In passing let us note that for a symmetric environment pinning can be proved, 
at least on the level of e(A) 0- Symmetric means that P^'^^ is concentrated on point 
configurations which are symmetric relative to the diagonal. In this case Ac = 1, i.e. 
e(A) = 2 for < A < 1 and e(A) > 2 for A > 1. Indeed in the symmetric case, the 
system is amenable to exact computations in terms of Fredholm determinants; a 
trace of the simplification can also be seen in our later discussion (see Section IT^ . 



2.2. Construction of the broken-lines process 

Hammersley gave a representation of the longest increasing subsequence problem 
for a random permutation in terms of broken lines built from a Poisson point process 
in the positive quadrant (we describe the construction in some detail below). The 
length of the longest increasing subsequence can then be seen to be the number of 
lines which separate the points (0, 0) and (n, n) from one another. The purpose of 
that representation is to obtain a superadditivity property which easily implies the 
existence of the limit in (|2.1|l — but doesn't specify its value. It is a very convenient 
formalism, which was used in 

[HElEi and [i3,[i3 

The broken line process Ts in a finite domain S can be defined as the space-time 
trace of some particle system with birth, death and immigration. For convenience 
we rotate the whole picture by an angle of tt/4 clockwise, so that the geodesic 
is restricted to never have a slope which is larger than 1 in absolute value. In 
what follows we will consistently use the letters t and x for the first and second 
coordinates, respectively, in the rotated picture; we will refer to t as "time" (the 
reason for that will become clear shortly) . The geodesies can then be seen as curves 
of space- type (using the usual language of general relativity). 
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Fig 1 . The construction of a Hammersley process. 

Let S be the planar, bounded domain defined (cf. Figure ^ as 

S ~ {{t,x) : to<t<h, g^{t) <x< g+it)}, (2.3) 

where to < ti are given points and g-(t) < g+{t), to < t < ti are piecewise hnear 
continuous functions such that, for some Iq-^ (resp. t^^) in the interval [to, ti], g+{t) 
(respectively g^{t)) increases (respectively, decreases) on (tg, tj^) (resp. (to, ^oi)) 
and decreases (respectively, increases) on {tol^ ti) (resp. {tol^ ^i)), always forming 
an angle of ±7r/4 with the t-axis. 

Consider four independent Poisson processes no,+ , Hq,-, II^. and n_ on the 
boundary of S. The processes no,± are supported on the leftmost vertical boundary 
component AqS" := {{to,x) : go,~ < x < go,+} of S, where go,± 9±{to), and 
they both have intensity ■\/ A'^^ /2. The process 11+ is defined on the "northwest" 
boundary A+5 := {{t,go,+{t)) : to < t < t^-^}, and the process n_ is defined on 
the "southwest" boundary A_5 := {{t,go,~{t)) : to < t < t^^}; their intensities 
(with respect to the length element of A ±5) are both V>J^. Finally, let Hin be 
a Poisson point process in S, with intensity X^^\ and independent of the previous 
four. 

Following the general strategy for the definition of Markov polygonal fields of 
0, we define a broken line process as follows. Each point of the Poisson process Hin 
is the point of birth of two particles which start moving in opposite directions, i.e. 
with velocities +1, —1. At each random point of no,+ , n„ a particle is born having 
velocity +1. Similarly at each random point of Ho.-, n^, a particle is born having 
velocity —1. All particles move with constant velocity until two of them collide, 
after which both colliding particles are annihilated (see Figure 

The state space Xs of the process Fs is the set of aU finite collections (71 , . . . , 7^) 
(including the empty one) of disjoint "broken lines" 7^ inside S. By a broken line in 
S we mean the graph "f = {{t, x) £ S : t = 7(x)} in 5* of a continuous and piecewise 
linear function 7, with slopes all ±1. Let Vs denote the probability distribution on 
Xs corresponding to the broken line process Ts defined above. 
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Let S" C S" be two bounded domains of the form of (|2.3|l . and let Vs' and Vs" be 
the probabiHty distributions of the broken line processes in S' and S" , respectively. 
The probability measure Vs" on Xs" induces a probability measure Vs" \s' on Xs' , 
which is the distribution of the restricted process Ts" H S". Then, by the choice 
we made of the boundary conditions, the following consistency property holds (see 

0): 

Vs"\s'=rs'. (2.4) 

This guarantees the existence of the broken line process F = T]^2 on R^, its dis- 
tribution V on X^2 being such that for every S of the above shape. Vis = Vs- 
Moreover V is invariant with respect to the translations of R^. 

Remark. The same description holds for the polynuclear growth (PNG) model, 
which describes a crystal growing layer by layer on a one-dimensional substrate 
through the random deposition of particles. They nucleate on existing plateaus of 
the crystal, forming new islands. In an idealization these islands spread laterally 
with constant speed by steady condensation of further material at the edges of the 
islands. Adjacent island of the same level coalesce upon meeting and on the top of 
the new levels further islands emerge. 

Observe that a path tt which can move only in the northeast-southeast cone (i.e. 
a path of space type) can collect at most one initial Poisson point from each broken 
line. In other words, the broken lines "factorize" the points of the configuration oj, 
in such a way that it tells us which points cannot be collected by the same path, 
and that the maximal number of points is bounded by the number of broken lines 
which lie in between start-point and end-point of the path tt. 

In fact, the lines also provide an explicit construction of a geodesic, as follows: 
Start at point {ti, x) on the right boundary and move leftward until you meet a line, 
then follow this line until you arrive at a point, which you can collect. Then start 
moving leftward again until you collect a point in the second-to-last broken line, 
and so on. The number of collected points is then essentially equal to the number 
of broken lines, though a little care needs to be taken as far as boundary conditions 
are concerned if this comparison is to be made completely formal. This observation 
led to a new proof of H2.1|l in |[||. It is also the starting point of our argument. 

2.3. Essential and non-essential points 

We now return to the question asked at the beginning of this section: How are extra 
added points affecting the initial system? We begin with a few purely determin- 
istic observations and statements. We will need some extra notations: Given any 
configuration w„ (not necessarily sampled from a Poisson process) of points in Q„, 
let H{ujn) be the number of broken lines produced by the above construction. For 
X = {x,t) and A = {xi, . . . , x^} we will denote by tl;„ U x or a;„ U A the configu- 
ration obtained from cj„ by the addition of the points xi, . . . , x^, and by F„(a;„) 
the associated configuration of broken lines. 

Proposition 2.1 (Abelian property, see For any choice oftUn, Xi and 

X2 we have 

F„ (i;„ U (xi U X2)) = F„ ((5„ U xi) U xj) = F„ ((ii„ U Xa) U Xi) . (2.5) 

Proposition 2.2 (Monotonicity, see jl7| ^. For any choice of ujn and xi we 

have 

i?(w„) < iJ(5„ U xi) < iJ(5„) + 1. (2.6) 
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Definition 2.3. Given ujn and x we say that x is essential for tj„ if H{LUn U xi) = 

Remark. The above definition is domain-dependent: If 5„ C Sm arc two domains 
in , and w„ , ojm are the restrictions of a configuration to to Sn and Sm respectively, 
then an extra added point which is essential for a;„ might not be essential for Um- 

If an added point x is essential, its presence is felt on the boundary of the domain 
by the appearance of an extra broken line going outside of the area. If an added 
point X is not essential, its presence can be felt on the boundary or not, but in any 
case it will change the local geometry of existing broken lines. 

Speaking informally, the configuration u5„ determines a partition of the domain 
into two (possibly disconnected) regions E and B, such that any additional point 
chosen in E will be essential for a;„, while if it is in B it will be not essential for 
uJn- It is easy to construct examples of configurations to for which E is empty (i.e., 
that are very insensitive to local changes). On the other hand B is never empty as 
soon as uj is not empty. It is also easy to give examples of the following situations: 

• xi is not essential for a;„ and X2 is not essential for cj„, but xi is essential for 
23,1 U X2 and X2 is essential for tZ;„ U Xi ; 

• xi is essential for a;„ and X2 is essential for ujn, but X2 is not essential for 
w„ U xi . 

Remark. In order to simplify our explanations and make some concepts as well as 
the notations more transparent (and lighter) , we will change from the consideration 
of point-to-point case to the point-to- hyperplane case, i.e. instead of looking for a 
geodesic connecting (0, 0) to (t, 0), we will be looking for an optimal path connecting 
(t, 0) to the line x = 0. This change is of purely "pedagogical" nature: All the 
ideas discussed above and below are easily transferred to the point-to-point case. 
Nevertheless we will not deny that it requires some amount of additional work due 
to boundary conditions. 

The reader should also not be surprised with our taking of starting point as (i, 0) 
and moving backward to the cc-axis in the point-to-hyperplane case (see the remark 
at the end of Section I2.2|l . Since our broken lines were constructed by drawing the 
space-time trajectories of particles in "forward time" , the information provided by 
the broken lines about the underlying point configuration is useful in the backward 
direction, and thus it forces us to construct the geodesic this way. Conversely, in 
order to construct a forward geodesic, one could construct the broken lines on the 
same point configuration but backward in time. 

Due to that, the area to which we will be restricting our process will be the 
triangular area Sn enclosed by segments connecting points (0, —n), (0, n) and (n, 0). 
If confusion doesn't arise, we will keep denoting related quantities by the same 
sub-index n as for the square case, for example uJn will stay from now on for the 
configuration of points in this triangular area. 

The next proposition is the crucial point in our construction. 

Proposition 2.4 (see |l7|). // x is essential for uJn then the point-to-plane 
geodesic in configuration a;„ U x has to collect point x. 

Again this is a purely deterministic statement, and does not depend on the choice 
of Zjn or x. Further we will be considering only cases when extra points are added 
only along the t-axis. 
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2.4- Propagation of influence 

Once an extra point x is added to the system, we need to update the configuration 
of broken fines. One way to do that is to redo the whole construction from scratch, 
i.e. to erase ah the existing broken fines and redraw them using the algorithm we 
described previously, taking the new point into account. It is then natural to ask 
how much the new picture differs from the old one, which is not perturbed by 
addition of point x. It turns out that there is a very simple algorithm allowing us 
to trace all the places of the domain where the addition of x will be felt, i.e. where 
local modification will be done. 

Consider an augmented configuration = a;„ U {x}. In order to see where 
and how the broken lines of r„(cj„) will be modified, we look at a new interacting 
particle system, starting from the points of uj'„, but with new interaction rules: 

1. Particles starting from the points of w„ are following the same rules as before, 
i.e. they are annihilated at the first time when they collide with any other 
particle. These particles will be called "regular particles" ; 

2. The two particles which start from the newly added point x (also with veloc- 
ities +1 and —1) will be called "superior particles", and they obey different 
rules: 

(a) Superior particles annihilate if and only if they collide with each other; 

(b) If a superior particle collides with some regular particle, the velocity of 
the superior particle changes to that of the incoming regular particle, 
which is annihilated while the superior particle continues to move. 

We will denote by p+ (resp. p~) the superior particle which starts from x with 
initial velocity +1 (resp. —1); The space-time trajectories of p+ and will be 
denoted by tt^ and tt^ , respectively. 

Observe that if any of these two trajectories leaves the triangular area Sn, it 
will never come back to it. (Notice that since superior particles can change their 
velocities during their evolution, both particles can leave the triangular area from 
the same side.) If 7r+ and intersect inside of Sn, then, according to rule I2al the 
corresponding superior particles are annihilated. 

The path 7r+ (resp. tt~) can be represented as an alternating sequence of con- 
catenated segments 7r+(l,-|-), 7r+(l,— ), 7r+(2,-|-), 7r+(2,— ),... (resp. 7r~(l,— ), 
TT^{1, +), (2, — ), (2, +), . . . ), where each segment corresponds to the time in- 
terval between two consecutive velocity changes of the superior particle p+ (resp., 
and during which its velocity is equal to -1-1 or —1 according to the sign given 
as second argument in the notation. 

The trajectories of p^ and p~ are completely determined by a;„ U {x}. Observe 
that each time p^ changes its velocity from -f-l to —1 (resp. p~ changes its velocity 
from —1 to +1), it starts to move along a segment which also belongs to some 
broken line 7^ from r„(tD„), and when it changes its velocity back to -1-1, it leaves 
this broken line, and moves until the next velocity flip, which happens exactly when 
the superior particle collides with the next broken line "/i+i in r„(a;„). This gives 
an extremely simple rule how to transform r„(cj„) into Tn{u)n) (see Figure EJ: 

• Erase all the tt^ (j, t) {i.e., all parts of the path of the superior particle which 
are contained in one of the original broken lines) to obtain an intermediate 
picture r^(w„); 

• Add all the tt^ (j, ±) thus obtaining the picture r"(a;„). 
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Fig 2. Propagation of influence. 



It is not hard to conclude that r"(cj„) = r„(cj^). In other words the paths 7r+ 
and TTjjT show how the "influence" of x spreads along the configuration r„(cD„). If 
both superior particles ^^^^ Px collide inside the domain Sn, the trajectories tt^ 
and 7r~ close into a loop, and outside of this loop the configuration r„(a;„) was not 
modified, i.e. the presence of x was not felt at all. 

Definition 2.5. Given a configuration lo of the underlying Poisson process, and an 
added point x on the time axis, we will denote by T(x;a;) (or simply r(x) if there 
is no confusion possible) the self-annihilation time of the pair of particles created 
at X, i.e. the time at which the paths tt" and 7r+ meet, if such a time exists; let 
r(x) = +00 otherwise. 

In the specific case we are looking at, t(x) is almost surely finite if A'^^' > 0, but 
it need not be the case for other underlying point processes. 

2.5. Interaction and Attractors 

Another important step in the analysis of the spread of influence, is to understand 
how the influence paths interact with each other if we add multiple points Xi , . . . , x^ 
to the initial configuration. Proposition 12.11 implies that we can obtain the full 
picture by adding the points one by one; to simplify the notations, in our description 
of the procedure we will also use the fact that the additional points will be placed 
along i-axis, though this is not essential. 

Again, due to the presence of time orientation, the nature of the interaction 
between influence paths becomes exposed in a more transparent way if we pro- 
ceed backward, i.e. if wc begin to observe the modifications first when adding the 
rightmost point, and then continue progressively, adding the points one by one, 
moving leftward, each time checking the effect created by the newly added point. 
For notational convenience let us index the new points in the backward direction, 
i.e. Xi = {ti, 0) with ti > t2 > ■ ■ ■ . 

Applying the construction described in the previous subsection successively for 
each of the new points, we obtain the following rules for the updating a configuration 
with multiple points added: Take the initial configuration ojn and let ujn — a;„ U 
{xi, 1 < i < r} be the modified configuration. In order to see where and how the 
broken lines of r„(w„) will be updated, consider a new particle representation built 
on the configuration ujn , and obeying the following rules: 

1. Regular particles, starting from the points of a;„, are annihilated as soon as 
they collide with any other particle; 
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2. The particles starting from the {:x.i)i<i<r, with velocities +1 and —1 are de- 
noted by and p~. respectively; again we shall call them superior particles. 
They behave as follows: 

(a) Whenever two superior particle of different types collide ("H — " colli- 
sion), they annihilate and both disappear; 

(b) If two superior particles of the same type collide 1-" or " " colli- 
sion), then they exchange their velocities (elastic interaction) and con- 
tinue to move; 

(c) If a superior particle collides with a regular particle, the velocity of the 
superior particle changes that of the incoming regular particle; the reg- 
ular particle is annihilated, while the superior particle survives and con- 
tinues to move. 

Denote the space-time trajectories of superior particles pj. and p^. by 7r+. and 
respectively. As before, each pair of paths n^. and Tr~. can be represented as 
an alternating sequence of concatenated segments 7r+(l,+), 7r+(l,— ), 7r+(2,+), 
7^^.(2, -),■•■ or TT" (1, -), TT" (1,-h), TT" (2, -), tt" (2, +),..., respectively, with the 
same convention for the sign of the velocities. We are now ready to complete the 
set of rules which govern the transformation of r„ ((!;„) in to Tn{ujn^): 

• Erase each segment tt^. {j, =p) from r„(w„), producing an intermediate picture 

• Add the segments tt^. (j, ±) to r^(a;„), thus producing rj((w„) = Tniujn^)- 

(Here as previously, two ± in the same formula are taken to be equal signs, while 
± and =F in the same formula stand for opposite signs.) 

Recall that we are working in the bounded triangular domain Sn with the configu- 
ration 2„U{xi, 1 < ^ < r}, where r is the number of added points. By fj^ = {t^ , ) 

and {~ = {tf- , Xf- ) we shall denote the end-points of the influence paths 7r+ and tt~. 
— they can be points where the corresponding paths exit the triangular domain, 
or points where a "H — " collision happens, in which case the two corresponding 
end-points are equal. 

Besides, let ~ {t^+ , 0) and = {t^- , 0) be the points of first return to the 

<-axis of the paths 7r+. and tt^. , respectively, and define 

:= min{tj+ , t^- , ,t^-}. (2.7) 

Last, let := 7r+ fl {{ti,x), x & R} and e~ := 7r~ n {{ti,x), x S M}, and denote 
by TT^. , 7r~. the parts of 7r+. and 7r~. lying between x^ and , and respectively 
between x,; and e,^ . 

Definition 2.6. Let Ji be the (random) Jordan curve starting at x^, following the 
path TT ~. until the point e~ , then the vertical line t — ti up to , and then the 
path 7?+. backward until it comes back to x^. The domain bounded by Jj will be 
called the attractor of the point Xj and denoted by Ai (see Figure EJ. The part of 
its boundary which is contained in tt,^. (resp. tt~.) will be called the upper (resp. 
lower) boundary of the attractor. 

It is important to understand how attractors are affected by one another, in 
order to give a convenient description of the whole augmented process. The key 
remark is the following: Informally speaking, a superior particle of a given type is 
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Fig 3. Attractors (shaded). 



not affected by an older one (i.e. one with larger index) of the same type. Indeed 
this is a consequence of the previous construction, and especially of the Abelian 
property (see ProDOsition l2.1|l . 

Of course this does not mean that the attractors Aj for j < jo do not change 
when is added, since it remains possible that a "H — "-collision happens. In 
that case, we get some kind of a monotonicity property, for the statement of which 
some additional notation will be needed. Recall that t{x;uj) stands for the self- 
annihilation time of a particle in the underlying scenery, i.e. ignoring the effect of 
the other new particles, both younger and older. We will denote by T(x,y;a;) the 
annihilation time of the set of superior particles born from x and y, i.e. the last 
time at which any of the corresponding four superior particles is still alive. Such an 
annihilation can happen in one of four ways: 

• Flat: X is born, the two particles thus created annihilate, then y is born and 
its two particles collide; 

• Embedded: x is born, then y appears between p~ and then the particles 
issued from y collide, then so do those issued from x; 

• Parallel: x is born, then y appears outside of (Pxi-Px); then the particles 
issued from y collide, then so do those issued from x; 

• Crossed: x is born, then y (also outside of (p~,p+)), then one particle issued 
from x annihilates the particle of the other type coming from y, then the 
remaining two collide. 

The combinatorics become much more involved when more particles are added; 
nevertheless, it is possible (if a bit technical if a formal proof is needed, see 0) to 
show the following: 

Proposition 2.7 (Monotonicity of the influence). For any two added points 
X and y, we have the following inequality: 

r(x,y;tj) > Max(T(x, cj), T(y; w)); 

and more generally the annihilation time of the union of two finite families of added 
points is at least equal to the larger of the two annihilation times of the parts. 

In the flat, embedded and parallel cases, the monotonicity extends to the shapes 
of the attractors (the attractor of y in the presence of x contains the one without); 
there is true reinforcement in the embedded case, in that the inclusion is strict as 
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soon as there is a "++"- or " "-collision. This is not always the case in crossed 

configurations (cf. Fig. [SJ where the attractor of X3 is shortened by the addition of 
X4), which leads us to the following definition: 

Definition 2.8. We say that two attractors Ai and Aj,i > j are connected if there 
exists a sub-sequence j = iq < ii < 12 < ■ ■ ■ < ik = i such that x^^ G -^v+i for all 
< r < fc. We will call {ij)o<j<k a connecting subsequence between Ai and Aj. 

Observe that if j > j > fc, if Ai is connected to Aj and Aj is connected to 
Ak, then Ai is connected to Ak- Nevertheless, due to the presence of orientation in 
the temporal direction, the above implication generally does not hold without the 
condition i > j > fc. 

Our construction immediately implies the following: 

Proposition 2.9. If Ai is connected to Aj, i > j, then ti > tj. 

Corollary 2.10. If Ai is connected to Aj, i > j, and the end-points and ej 
belong respectively to the south-east and north-east boundaries of the triangular 
domain Sn, then so do and e^. 

Corollary 2.11. Assume that Ai is connected to Aj, i > j, with connecting sub- 
sequence {is)o<s<k-' If ^j is essential for the configuration a;„, then so are all the 
Xi, , 1 < s < fc. 

2. 6. Pinning of the geodesies 

We now return to our original problem. Observe that if, for some fixed configuration 
u)n in Sn, we pick a realization of the points (xi)i<i<t, in such a way that Ai and Av 
are connected (say), with connecting subsequence {ij)o<j<k, then all the x^^. must 
be essential for u;„, and therefore the point-to-plane geodesies for the configuration 
Un U {xi} has to visit all the x,;^ . 

In the new formalism, the original question of whether, for any given density 
A(i) > of the one-dimensional Poisson point process, the limiting value in ()2.1|l 
is increased, becomes equivalent to the following: Is there a positive S such that, 
with high probability as n goes to infinity, at least a fraction S of the newly added 
points are essential for w^^-'? 

This question is more complicated than simply whether there exists a chain of 
pairwise connected attractors spanning from the left to the right boundary of the 
domain: Indeed, such a chain does not necessarily have positive density. In the next 
Section we also mention some of the interesting mathematical questions that arise 
in the construction. 

It is not an easy task to understand how the attractors behave. The fact that the 
structure of the influence paths 7fxi+i and tt".^^ depends on r„(a;„U{xj, 1 < j < i}), 
but not on the x^ for j > z -I- 1, reduces the problem to checking whether none of 
the influence paths TTxi+iJ ^xi+i ^i^s the t-axis before x^, in which case Ai+i is 
connected to Ai (see Figure 

For a single point x added to the initial configuration, each infiucnce path 7r+, 
7r~ has the same statistical properties as what is known as a "second-class particle" 
in the framework of exclusion processes. Since in the definition of an attractor an 
important role is played by the (possible) return times of the infiucnce paths to the 
t-axis, several things must be settled: 

1. The first return time to the t-axis of a single influence path. It is believed (but 
remains a challenging open problem) that in the case of a one-dimensional 
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exclusion process, a second-class particle behaves super-diffusively. Though 
some bounds are available, and we know the mean deviations of the second 
class particle they do not provide good control on return times; 
2. The joint behavior of the influence paths tt^ and tt^. Generally, it is a com- 
plicated question too, but for our purposes wc need to have such control only 
up to the first times when one of tt^, ti~ returns to the i-axis. Before such a 
time, both paths stay apart from each other, and some good mixing proper- 
ties of the system come into play; so the question reduces to how efficiently 
we control point ^ 

The fact that the influence lines of "younger" points (with smaller indices, i.e. 
sitting more to the right) repeal the influence lines of "older" points, leads to 
the following observation: Once the attractor Ai+m of an older point reaches the 
younger point x^, then it cannot end before the attractor At ends. If the attractor 
Ai ends before reaching the next point Xi_i, then Ai+m can still go forward, and 
possibly itself reach Xi_i. Observe that at the time Ai ends, the boundaries of Ai+m 
are necessarily at a positive distance from the i-axis. 

That, together with the fact that the evolution of p^. and p~. in the slab 
(ti, ti-i) X R depends only on r„(il;„) brings some notion of week dependence to the 
system from one side, and the idea of a "re-start point" from another. This reduces 
the study of percolation of attractors to a more general problem of one-dimensional, 
long-range, dependent percolation which wc formulate in the next section. There 
we also mention some of the interesting mathematical questions that have arisen. 

3. Stick percolation 

In this section we introduce two "stick percolation" models, which will serve as toy 
models in the study of the propagation of influence in the broken-line model. In 
spite of their apparent simplicity, these models can be very useful studying effects of 
columnar defects and establishing (bounds for) critical values for asymptotic shape 
changes for some well known one-dimensional growth systems (see |^). 

3.1. Model 1: overlapping sticks 

Let {xi)i^fi be a Poisson point process of intensity A > on the positive real 
line. We call the points of this process "seeds" , and assume that they arc ordered, 
Xq being the point closest to the origin. To each seed Xi we associate a positive 
random variable Si (a "length") and assume that the {Si)i^iq are i.i.d. with common 
distribution function F. 

The system we consider is then the following: For every i G N, construct the 
segment Si = [xi,Xi + Si] (which we will call the i-th "stick"). We say that the 
sticks Si and Sj, i < j are connected if Xj < Xi + Si, i.e. if they have non-empty 
intersection; the system percolates if and only if there is an infinite chain of distinct, 
pairwise connected sticks, which (with probability 1) is equivalent to saying that 
the union of the sticks contains a half-line. 

It is easy to see that the system percolates with probability or 1 (it is a tail 
event for the obvious filtration); and in fact there is a complete characterization of 
both cases: 

Proposition 3.1 (|4!). Let R{x) = P{Si > x) be the tail of the stick length 
distribution, and let ip{x) — R{u)du. Then the system percolates with probability 
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1 if and only if 




e-^'^^^'dx < +00. 



(3.1) 



This leads us to the foUowing definition: 

Definition 3.2. The distribution F governing the stick process is said to be cluster- 
stable if the system percolates for every positive value of A. 

Example 1 (Return times of random walks). One natural distribution for 
the length of the sticks, in view of the previous construction, is the following: At 
time Xi, start two random walks with no drift (the specifics, e.g. whether they are 
discrete or continuous time walks, will not matter at this point — for that matter 
we could also take two Brownian motions), one from +1 and the other from — 1. 
Then, let Xi + Si be the first time when these two walks meet. It is well known that, 
up to multiplicative constants, P{Si > t) behaves as t~^^^ for large t, as soon as 
the walks are irreducible and their step distribution have finite variance. It is easy 
to check that the obtained distribution is cluster-stable. 

Example 2. Change the previous example a little, as follows: For every value of i, 
start a two-dimensional Brownian motion (or random walk) starting at (x^.i), and 
let Xi + Si be the first hitting point of the axis by this Brownian motion; but erase 
the stick if 5*^ < 0. Then the distribution F is Cauchy restricted to be positive, and 
its tail is equivalent to P{Si > t) ^ c/t as t goes to infinity. In that case, the stick 
length distribution is not cluster-stable for model 1, and there exists a critical value 
Ac = 1/c, such that for A > Ac the system percolates while it does not if A < Ac. 

3.2. Model 2: reinforced sticks 

There are several ways to mimic the interaction between "funnels" . One of the most 
simple possibilities is to add extra rules to Model 1 to account for the interaction. 
The basic idea behind this modification, is that if two or more sticks from Model 
1 overlap, then there is a certain reinforcement of the system, which depends on 
how many sticks overlap, and then the whole connected component is enlarged 
correspondingly. 

One way to do that can be described informally as the following dynamic process. 
First, see each stick Si as the flight time of a particle tt^ born at time Xi. We want to 
model the fact that if tt^ wants to land when a younger one (say ttj) is still flying, 
instead tt^ "bounces" on nj; ttj on the other hand should not be affected by tt^, 
if the process is to look like the propagation of influence described in he previous 
section. 

Assign to each particle a "counter of chances" Ni which is set to 1 when the 
particle is born (formally it should be a function from R+ to N, and we let Ni{xi) = 
1). Then, two things can happen. If Xi -\- Si < there is no interaction and tt^ 

dies at time Xi + Si. If on the other hand Xi + Si > at time .t,;+i the particle 

TTi gets a "bonus", so that Ni{xi+i) = 2; and similarly, it gets a bonus each time it 
passes above a seed point, so that Ni{xi -\-Si—) — l is the number of seed points 
in Si. Now when tt lands, its counter is decreased by 1, but if it is still positive the 
particle bounces on the axis and restarts using an independent copy of S. If all the 
chances are exhausted before the closest seed is reached, Tr^ gets killed. 

Again, we may ask a similar question: Given A > 0, for which distribution func- 
tions F is the probability for a given particle to survive up to infinity positive? It is 
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obvious that even with £'(5') < +00 we can obtain infinite trajectories for certain 
(large) values of A. We will call F cluster- stable for the reinforced process if this 
happens for every positive A. 

Example. The Cauchy distribution used in the previous example cluster-stable for 
the reinforced model. 

It is easy to see {e.g. by a coupling argument) that if the original system perco- 
lates, the reinforced version (for the same value of A and the same length distribu- 
tion F) percolates too. In particular a cluster-stable distribution is cluster-stable for 
the reinforced problem. Nevertheless it is an interesting open problem to give full 
characterization of distributions which are cluster-stable for the reinforced models. 
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